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High Fidelity Low-Order Modeling For Multi-Staged 
Separation Processes Based Upon The Compartmental Approach 

1 . Field of the Invention 

This invention relates to process systems engineering 
and more particularly to the simulation of separation 
processes in the process industries. 

2 . Description of the Prior Art 

Process engineering involves the design of a wide 
variety of processing plants and processes carried out 

10 therein. Such processes include, but are not limited to, 
chemical, petrochemical, refining, pharmaceutical, and 
polymer industries. In process engineering, the use of 
computer based models to develop and evaluate new 
processes, design and retrofit plants, and optimize the 
operation of existing plants is rapidly becoming a 
standard. At every stage of process design, development 
and operation, rigorous models generated by process 
simulation software systems can be used to make better 
engineering and business decisions. 

20 In a process simulation software system, the 

performance of a process industry plant in which there is 
a continuous flow of materials and energy through a 
network of process units is simulated. The process unit 
can include equipment, such as distillation columns, 
retaining vessels, heating units, pumps, conduits, etc. 
Typically, the process simulation software features 
computer models that allow process engineers to simulate 
the operation of various pieces of equipment used in a 
proposed or existing manufacturing process. The end 

30 results from the simulation software system provide a 
display of the simulated performance of the plant under 
various conditions and estimate of the capital and 
operating cost of the plant and its profitability. 

Multi-staged separation processes are used in the 
chemical industry in applications such as distillation 
columns. The prior art has not found it easy to solve 
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the dynamic response of such processes in a reasonably 
short time through modeling while maintaining model 
fidelity for on-line application purposes. 

In the most rigorous applications of dynamic 
simulation of stage-wise separation processes, high order 
models are usually employed. These models typically 
involve the solutions of material and energy balance 
equations for every stage of the process. The high 
dimensionality of these models, however, make them 
10 impractical in those applications such as on-line 
optimization and control, and dynamic data reconciliation 
where time consideration is important. In these 
applications, solutions must be obtained within a limited 
time frame. 

Several low order methods that involve significantly 
fewer differential equations, such as linearization with 
model reduction and orthogonal collocation, have been 
proposed in the engineering literature. However, they are 
all deficient in some respects. For example, 

20 linearization with model reduction would render the model 
parameters with no physical significance, while the 
orthogonal collocation approach would give results that 
may not match the discrete nature of real world problems. 

The academic research project reported by Benallou 
A., Seborg, D.E., and Mellichamp, D.A., "Dynamic 
Compartmental Models For Separation Processes", AIChE 
Journal, vol 32(7), pp. 1067-1078 (July, 1986) resulted 
in the development of the compartmental approach for 
distillation column applications. The academic research 

30 projects reported by Horton, R.R., Bequette, B.W., and 
Edgar, T.F., "Improvements In Dynamic Compartmental 
Modeling For Distillation", Computers and Chemical 
Engineering, vol. 15(3), pp. 197-201 (1991); and Musch, 
H.E. and Steiner, M . , "Order Reduction Of Rigorous 
Dynamic Models For Distillation Columns", Computers and 
Chemical Engineering, vol. 17S, pp. S311-S316 (1993), 
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have resulted in refinements to the compartmental 
approach. All three of these methods, however, suffer 
from the following deficiencies: 

1. Difficulties in initialization; and 

2. Difficulties in converging to a solution. 

The present invention addresses these deficiencies 
and thus adapts more readily to real world problems. 

Summary of the Invention 
A medium having instructions for performing steps for 
10 creating a regular model for a multi-stage separation 
process having feed trays, reboiler return trays and 
condenser return trays and one or more additional trays 
that allows the regular model to converge to a solution 
comprising: 

(a) lumping each of the one or more additional 
trays into an associated compartment that comprises one 
or more of the one or more additional trays that are not 
in another associated compartment; 

(b) treating the feed trays, reboiler return trays 
20 and the condenser return trays as one tray compartments; 

(c) designating one of the one or more additional 
trays in each of the associated compartments comprising 
one or more additional trays as a sensitive tray; 

(d) equating liquid holdup in each of the 
compartments to the total liquid holdup of the trays in 
each of the compartments; 

(e) ignoring vapor holdup in each of the 
compartments; 

(f) calculating phase equilibria in an associated 
30 compartment only for the sensitive tray; 

(g) presuming that all trays other than the 
sensitive trays in each of the associated compartments 
respond instantaneously to liquid and vapor flow; and 

(h) basing temperatures of all trays other than the 
sensitive trays in each of the associated compartments on 
linear interpolation between the temperature at each of 
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two adjacent sensitive trays. 

A medium having instructions for performing steps 
for creating an initialization model for a multi-stage 
separation process having feed trays for receiving a feed 
stream flow, reboiler return trays for providing a liquid 
flow to a reboiler and receiving from the reboiler a 
vapor flow and condenser return trays for providing a 
vapor flow to a condenser and receiving from the 
condenser a liquid flow and one or more compartments with 
multiple trays having vapor and liquid streams flowing 
therethrough that generates a set of initial values for a 
regular model of the process, comprising: 

(a) converting all of the differential equations in 
the regular model to steady state equations by setting 
the derivative term of all of the differential equations 
to zero; 

(b) taking the temperature of the condenser return 
tray as the average of the temperatures of the liquid 
flow from the condenser and the feed stream flow; 

(c) taking the temperature of the reboiler return 
tray as the average of the temperatures of the vapor flow 
from the reboiler and the feed stream flow; 

(d) taking the liquid flow in the condenser return 
tray to be the same as the liquid flow from the 
condenser; 

(e) taking the vapor flow in the reboiler return 
tray to be the same as the vapor flow from the reboiler; 

(f) taking the liquid flow in the feed tray as the 
sum of the liquid flow from the condenser and the feed 
stream flow; 

(g) assuming that all of the trays have 100% 
efficiency; and 

(h) assuming there is no interaction between the 
vapor and liquid flows through the compartments with 
multiple trays. 

A method for creating a regular model for a multi- 
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stage separation process having feed trays, reboiler 
return trays and condenser return trays and one or more 
additional trays that allows the regular model to converge 
to a solution. The method is: 

(a) lumping each of the one or more additional 
trays into an associated compartment that comprises one 
or more of the one or more additional trays that are not 
in another associated compartment; 

(b) treating the feed trays, reboiler return trays 
and the condenser return trays as one tray compartments; 

(c) designating one of the one or more additional 
trays in each of the associated compartments comprising 
one or more additional trays as a sensitive tray; 

(d) equating liquid holdup in each of the 
compartments to the total liquid holdup of the trays in 
each of the compartments; 

(e) ignoring vapor holdup in each of the 
compartments; 

(f) calculating phase equilibria in an. associated 
compartment only for the sensitive tray; 

(g) presuming that all trays other than the 
sensitive trays in each of the associated compartments 
respond instantaneously to liquid and vapor flow; and 

(h) basing temperatures of all trays other than the 
sensitive trays in each of the associated compartments on 
linear interpolation between the temperature at each of 
two adjacent sensitive trays. 

A method of process engineering which comprises 
creating a regular model for a multi-stage separation 
process having feed trays, reboiler return trays and 
condenser return trays and one or more additional trays 
that allows the regular model to converge to a solution 
according to prescribed guidelines, 

the guidelines comprise: 

(a) lumping each of the one or more additional 
trays into an associated compartment that comprises one 
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or more of the one or more additional trays that are not 
in another associated compartment; 

(b) treating the feed trays, reboiler return trays 
and the • condenser return trays as one tray compartments; 

(c) designating one of the one or more additional 
trays in each of the associated compartments comprising 
one or more additional trays as a sensitive tray; 

(d) equating liquid holdup in each of the 
compartments to the total liquid holdup of the trays in 

10 each of the compartments; 

(e) ignoring vapor holdup in each of the 
compartments; 

(f) calculating phase equilibria in an associated 
compartment only for the sensitive tray; 

(g) presuming that all trays other than the 
sensitive trays in each of the associated compartments 
respond instantaneously to liquid and vapor flow; and 

(h) basing temperatures of all trays other than the 
sensitive trays in each of the associated compartments on 

20 linear interpolation between the temperature at each of 
two adjacent sensitive trays. 

A method for creating an initialization model for a 
multi-stage separation process having feed trays for 
receiving a feed stream flow, reboiler return trays for 
providing a liquid flow to a reboiler and receiving from 
the reboiler a vapor flow and condenser return trays for 
providing a vapor flow to a condenser and receiving from 
the condenser a liquid flow and one or more compartments 
with multiple trays having vapor and liquid streams 

30 flowing therethrough that generates a set of initial 
values for a regular model of the process. The method 
is : 

(a) converting all of the differential equations in 
the regular model to steady state equations by setting 
the derivative term of all of the differential equations 
to zero; 
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(b) taking the temperature of the condenser return 
tray as the average of the temperatures of the liquid 
flow from the condenser and the feed stream flow; 

(c) taking the temperature of the reboiler return 
tray as the average of the temperatures of the vapor flow 
from the reboiler and the feed stream flow; 

(d) taking the liquid flow in the condenser return 
tray to be the same as the liquid flow from the 
condenser; 

10 (e) taking the vapor flow in the reboiler return 

tray to be the same as the vapor flow from the reboiler; 

(f) taking the liquid flow in the feed tray as the 
sum of the liquid flow from the condenser and the feed 
stream flow; 

(g) assuming that all of the trays have 100% 
efficiency; and 

(h) assuming there is no interaction between the 
vapor and liquid flows through the compartments with 
multiple trays. 

20 A method of process engineering which comprises 

creating an initialization model for a multi-stage 
separation process having feed trays for receiving a feed 
stream flow, reboiler return trays for providing a liquid 
flow to a reboiler and receiving from the reboiler a vapor 
flow and condenser return trays for providing a vapor flow 
to a condenser and receiving from the condenser a liquid 
flow and one or more compartments with multiple trays 
having vapor and liquid streams flowing therethrough that 
generates a set of initial values for a regular model of 

30 the process, according to prescribed guidelines, 
the guidelines comprise: 

(a) converting all of the differential equations in 
the regular model to steady state equations by setting 
the derivative term of all of the differential equations 
to zero; 
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(b) taking the temperature of the condenser return 
tray as the average of the temperatures of the liquid 
flow from the condenser and the feed stream flow; 

(c) taking the temperature of the reboiler return 
tray as the average of the temperatures of the vapor flow 
from the reboiler and the feed stream flow; 

(d) taking the liquid flow in the condenser return 
tray to be the same as the liquid flow from the 
condenser; 

10 (e) taking the vapor flow in the reboiler return 

tray to be the same as the vapor flow from the reboiler; 

(f) taking the liquid flow in the feed tray as the 
sum of the liquid flow from the condenser and the feed 
stream flow; 

(g) assuming that all of the trays have 100% 

efficiency; and 

(h) assuming there is no interaction between the 
vapor and liquid flows through the compartments with 
multiple trays. 

20 A computer executable model for a multi-stage 

separation process having feed trays, reboiler return 
trays and condenser return trays and one or more 
additional trays. The model comprises: 

(a) each of the one or more additional trays lumped 
into an associated compartment that comprises one or more 
of the one or more additional trays that are not in 
another associated compartment; 

(b) the feed trays, reboiler return trays and the 
condenser return trays each treated as one tray 

30 compartments; 

(c) one of the one or more additional trays in each 
of the associated compartments comprising one or more 
additional trays designated as a sensitive tray; 

(d) liquid holdup in each of the compartments 
equated to the total liquid holdup of the trays in each 
of the compartments; 
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(e) vapor holdup in each of the compartments 
ignored; 

(f) phase equilibria in an associated compartment 
calculated only for the sensitive tray; 

(g) all trays other than the sensitive trays in 
each of the associated compartments are presumed to 
respond instantaneously to liquid and vapor flow; and 

(h) temperatures of all trays other than the 
sensitive trays in each of the associated compartments 
are based on linear interpolation between the temperature 
at each of two adjacent sensitive trays. 

A computer executable model having differential 
equations for a multi-stage separation process having feed 
trays for receiving a feed stream flow, reboiler return 
trays for providing a liquid flow to a reboiler and 
receiving from the reboiler a vapor flow and condenser 
return trays for providing a vapor flow to a condenser and 
receiving from the condenser a liquid flow and one or more 
compartments with multiple trays having vapor and liquid 
streams flowing therethrough. The model comprises: 

(a) all of the differential equations converted to 
steady state equations by setting the derivative term of 
all of the differential equations to zero; 

(b) the temperature of the condenser return tray is 
taken as the average of the temperatures of the liquid 
flow from the condenser and the feed stream flow; 

(c) the temperature of the reboiler return tray is 
taken as the average of the temperatures of the vapor flow 
from the reboiler and the feed stream flow; 

(d) taking the liquid flow in the condenser return 
tray to be the same as the liquid flow from the condenser; 

(e) taking the vapor flow in the reboiler return 
tray to be the same as the vapor flow from the reboiler; 

(f) taking the liquid flow in the feed tray as the 
sum of the liquid flow from the condenser and the feed 
stream flow; 
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(g) assuming that all of the trays have 100% 
efficiency; and 

(h) assuming there is no interaction between the 
vapor and liquid flows through the compartments with 
multiple trays. 

Description of the Drawing 
Fig. 1, the only drawing figure, shows a distillation 
column as one example of a multi-staged separation process 
to which the present invention can be applied. 

10 Description of the Preferred Embodiment (s) 

Fig. 1 shows a section of a distillation column 10, a 
typical staged separation process, in which a stream of 
fresh feed 12a is introduced into the column 10 and 
separated into groups of components. The column 10 
includes a number of horizontal trays 30a to 30o which 
provide the surfaces on which mass transfer (component 
exchange) can take place between the cold liquid 12b 
returning from a condenser (not shown) and flowing 
downward from the top of column 10 and the hot vapor 12c 

20 returning from a reboiler (not shown) and flowing upward 
from the bottom of column 10. 

As is shown in Fig. 1 and in accordance with the 
present invention, several trays such as those in 
compartments 14, 16, 18 are treated as a single entity for 
simulation purposes. Compartment 14 has trays 30b, 30c, 
30d, 30e; compartment 16 has trays 30f, 30g, 30h, 30i; and 
compartment 18 has trays 30k, 301, 30m, 30n. 

Column 10 also includes a top compartment 20 that has 
a single condenser return tray 30a. Compartment 20 

30 provides the hot vapor 12d flowing upward to the condenser 
and receives the cold liquid 12b returning from the 
condenser. Column 10 also has a feed compartment 22 that 
has a single feed tray 30 j which receives the stream 12a 
of fresh feed. Column 10 further has a bottom compartment 
24 that has a single reboiler return tray 30o. 
Compartment 24 receives the upwardly flowing hot vapor 12c 
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from the reboiler and provides the downwardly flowing cold 
liquid 12e to the reboiler. 

As is well known the distillation column 10 includes 
a rectifying section and a stripping section. The 
rectifying section consists of the trays directly above 
the feed compartment 22 and thus includes compartments 14 
and 16. The stripping section consists of the trays 
directly below the feed compartment 22 and thus includes 
compartment 18. 

10 The present invention provides guidelines for 

creating a regular model for a staged separation process 
which in this example is a distillation column. These 
guidelines should be used in the context of the basic 
separation model formulation known as the MESH equations. 
As is well known, the MESH equations are the mass and 
component balance equations, phase equilibrium equations, 
summation equations, and enthalpy (energy) balance 
equations. For rigorous modeling, the MESH equations are 
applied to each tray or section of a staged process. The 
20 guidelines of the present invention for creating the 
regular model for the process which are given below 
represent a deviation from the normal practice. The 
result of using these guidelines is a computer executable 
regular m odel of the multi-stage separation process that 
has all of the features called for in the guidelines. 

The guidelines for creating the regular model for the 
staged separation process are: 

1. Several trays are, as is described above, lumped 
into a compartment as is shown in Fig. 1 for 
30 compartments 14, 16 and 18. One of the trays in each 

of compartments 14, 16 and 18 is taken to be the 
"sensitive" tray of that compartment. The sensitive 
tray may be any tray in the compartment, for 
example, the middle tray in each of compartments 14, 
16 and 18. 
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2. Feed trays such as tray 30 j, reboiler return trays 
such as tray 30o, and condenser return trays such as 
tray 30a are not lumped. They are treated as one- 
tray compartments as is shown in Fig. 1 for 
compartments 20, 22 and 24. 

3. The liquid holdup in a compartment is equal to the 
total holdup of the constituent trays. As is well 
known the liquid holdup is the liquid that is 
consistently moving through a tray. 

10 4 . Vapor holdup is ignored. As is well known the 

response to a change in vapor is much faster than a 
response to a change in liquid. Tray delta pressure 
is a parameter to be estimated, for example, from 
overhead and bottom pressure measurements. 
5. Phase equilibria in a multi tray compartment are 
calculated only for the sensitive trays and not for 
all of the individual trays. As is well known phase 
equilibria depends on the pressure, temperature and 
composition of a mixture. 

20 6. The compositions, vapor flowrates, and liquid 

flowrates of the non-sensitive trays, that is the 
the trays in a multi-tray compartment other than a 
sensitive tray are based on steady state balances 
with respect to the sensitive tray, and with vapor- 
liquid equilibrium ratios the same as those for the 
sensitive tray. In other words, the non-sensitive 
trays are presumed to respond instantaneously to 
changes in liquid and vapor flow. 
7 . Temperatures of the non-sensitive trays are based on 

30 linear interpolation between the temperature at each 

of the two adjacent sensitive trays. 

The guidelines described above for creating the 
regular model for the process solves the problem of 
failure to converge to a solution that exists in the 
methods of the prior art compartmental approach to a 
distillation column. 
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To solve the problem of initialization failure, a 
separate computer executable model is used to generate a 
set of initial values for the regular model for the 
staged separation process. This computer executable 
initialization model is similar to the regular model as 
described above but with the following additional 
features : 

1. A11 differential equations that are in the regular 
model are converted to steady state equations by 

10 setting the derivative terms to zero. 

2. The top tray 30a temperature is taken to be the 
average of the condenser return and the feed 
temperatures. Likewise, the bottom tray 30o 
temperature is taken to be the average of the 
reboiler return and the feed temperatures. 

3. The top tray liquid flow 30a is taken to be the same 
as the condenser return flow 12d. The feed tray 30j 
liquid flow is the sum of the top and feed flows. 
The bottom tray 30o vapor flow is the same as the 

20 reboiler return flow 12e. 

4. Tray efficiency is assumed to be 100% for all trays. 

5. There is no interaction between the vapor and liquid 
streams in both the rectifying and the stripping 
sections . 

The guidelines of the present invention described 
above for setting up the regular model for the 
distillation column 10 and for setting up the 
initialization model for column 10 can also be applied to 
the modeling and simulation of other types of multi-staged 
30 separation processes, e.g., a packed tower. As is well 
known, a packed tower is a cylindrical column filled with 
a mass of inert solid shapes called tower packing. The 
liquid, distributed from above, uniformly wets the 
surfaces of the packing as it flows down the column. The 
vapor, fed from below, flows upward through the 
interstices in the packing. The packing thus provides a 
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large area of surface on which mass transfer can take 
place between the liquid and the vapor. 

A packed tower may be simulated as a distillation 
column with multiple trays. The equivalent number of trays 
in a packed column may be calculated from the column 
height and a parameter known as the height equivalent to a 
theoretical plate (HETP) . The HETP is usually in the 
range of 1 to 2 feet depending on the type of packing 
used. Therefore, a 30-ft packed column may be simulated 
10 as a distillation column with 15 to 30 trays. Once the 
number of equivalent trays is determined, the guidelines 
as described in the previous sections may be applied to 
the simulation of the packed column. 

Thus as can be appreciated from the description above 
the guidelines of the present invention can be used with 
any multi-stage separation process 

As can be appreciated the present invention, 
contrary to the prior art which stipulates that dynamic 
component balances are to be performed for all the trays, 
20 calls for steady-state component balances to be performed 
for the non-sensitive trays only. Further, instead of 
constant volatilities, constant vapor-liquid equilibrium 
ratios are used for the non-sensitive trays. These lead 
to faster and more robust numerical solutions for on-line 
applications . 

As can also be appreciated the combined effect of 
the features of the initialization model developed in 
accordance with the guidelines described herein allows 
the initialization model to generate a set of initial 
30 values for the regular model without which the regular 
model would often fail to initialize. 

It is to be understood that the description of the 
preferred embodiment (s) is (are) intended to be only 
illustrative, rather than exhaustive, of the present 
invention. Those of ordinary skill will be able to make 
certain additions, deletions, and/or modifications to the 



E20020390 ]_5 

embodiment (s) of the disclosed subject matter without 
departing from the spirit of the invention or its scope, 
as defined by the appended claims. 



